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Abstract 

We present a concise review of the recent important experimental developments on neutrino mixing (hints for sterile neutrinos, 
large O13, possible non maximal O23, approaching sensitivity on Sep) and their implications on models of neutrino mixing. The new 
data disfavour many models but the surviving ones still span a wide range going from Anarchy (no structure, no symmetry in the 
lepton sector) to a maximum of symmetry, as for the models based on discrete non-abelian flavour groups. 
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1. Introduction 









Quantity 


ain 2013 


sin 


T2K[6] 
MINOS [7] 

DC[8] 
RENO[9] 
DYB[10] 


ir° " 

041+0047 

0.0109 ±0.039 
0.113 ±0.023 
0.0890 ±0.01 12 


0.028!«:«^: 

0.028 ±0.011 
0.029 ± 0.006 
0.023 ± 0.003 



On the experimental side the main recent develop- 
ments on neutrino mixing [ 1, 2, 3, .4, were the results 
on 6»i3 from T2K|i], MINOsH, DOUBLE CHOOZjl, 
RENO and especially Day a Bay fl^ (see Table 1, 
where the latest numbers presented at the summer '12 
conferences were collected). The combined value of 
is by now about 10 cr away from zero and the central 
value is rather large. In turn a sizable On allows to ex- 
tract an estimate of 623 from accelerator data like T2K 
and MINOS. There are solid indications of a deviation 
of 023 from the maximal value, in the first octant. In ad- 
dition, some tenuous hints that cos 6cp < are starting 
to appear in the data. A summary of recent global fits to 
the data on oscillation parameters is presented in Table 

2 im, id, iH. 

On the possible existence of sterile neutrinos a num- 
ber of hints have been recently reported. They do not 
make yet an evidence but certainly pose an experimen- 
tal problem that needs clarification. First, there is the 



Table 1. The reactor angle measurements from the recent experi- 
ments T2K|6], MINOSI7], DOUBLE CHOOZ| 8], DayaBay 1 10] and 
RENO f9*], for the normal hierarchy case (in the inverse hierarchy case 
the values do not differ by much) 



MiniBooNE experiment that has recently published 111 
a combined analysis of v,, appearance in a beam to- 
gether with Ve appearance in a beam. An excess 
was observed for events with (anti)neutrino energy in 
the range 200 < Ey < 3000 MeV (3.8cr). The al- 
lowed regions resulting from a two-neutrino fit of the 
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data are consistent with v.. 



and v„ — > Vf oscilla- 
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Quantity 


Ref. fll] 


Ref. [12] 


Ami, (10-^ eV^) 
AmL (10-3 eV^) 

sin^ 012 

sin^ 023 

sin^ 013 


7 ^4+0.26 
o 4^+0.06 

o.307t;- 

0.386^0^4 
0.0241 + 0.025 


7.50 ±0.185 

7 47+0.069 
^•^ -0.067 

0.30 ±0.013 

41+0.037 
"■^^-0.025 

0.023 ± 0.0023 



Table 2. Fits to neutrino oscillation data. For sin 623 only the absolute 
minimum in the first octant is shown 



tions in the 0.01 to 1 eV^ range for Am^ and consistent 
with the allowed region reported by the LSND experi- 
ment. In both channels, but especially in the neutrino 
beam case, most of the excess is at low Ey, below 500 
MeV, where the evaluation of the background is par- 
ticularly difficult. Recently the ICARUS experiment at 
Gran Sasso has published the results of a search for elec- 
trons produced by the CERN beam [ 15]. No excess over 
the background was observed. As a consequence most 
of the region allowed by LSND, MiniBooNE. KAR- 
MEN... is excluded and only a small domain around 
(Am'-, sin2(20)) = (0.5ey^ 0.05)) is left. Then there is 
the reactor anomaly: a reevaluation of the reactor flux 
lfl6l l produced an apparent gap between the theoretical 
expectations and the data taken at small distances from 
reactors (< 100 m). The discrepancy is of the same order 
of the quoted systematic error whose estimate, detailed 
in the paper, should perhaps be reconsidered. Actually, 
a different analysis confirmed the normalization shift 
Iit Ii. Similarly the Gallium anomaly |18] depends on 
the assumed cross-section which could be questioned. 
Cosmological data allow, actually even favour, the ex- 
istence of one sterile neutrino, while the most stringent 
bounds arising form nucleosynthesis disfavour two or 
more sterile neutrinos (assuming that they are thermal- 
ized) Over all, only a small leakage from ac- 

tive to sterile neutrinos is allowed by present neutrino 
oscillation data (see, for example, refs. 11211 l22l 12311 
and references therein). If all the indications listed 
above were confirmed (it looks unlikely) then adding 
one or more sterile neutrinos would probably not be 
enough to satisfactorily describe all the data. There is 
in fact a strong tension between appearance (LSND, 
MiniBooNE, ICARUS) and disappearance (reactors, Ga 
anomaly) data. Thus the situation is at present very 
confuse but the experimental effort should be continued 
because establishing the existence of sterile neutrinos 
would be a great discovery. In fact a sterile neutrino 
is an exotic particle not predicted by the most popular 
models of new physics. In the following we restrict our 



discussion to 3-neutrino models. 

The rather large measured value of 0i3, close to the 
old CHOOZ bound and to the Cabibbo angle, and the 
indication that 023 is not maximal both go i n the di- 
rection of models based on Anarchy 124 12511 . i.e. the 
ansatz that perhaps no symmetry is needed in the lep- 
tonic sector, only chance (this possibility has been re- 
cently reiterated, for example, in Ref. ll26ll ). Anar- 
chy can be formulated in a 5 U{5) ® U(l)fff context by 
taking different Froggatt-Nielsen ll27ll charges only for 
the S U(5) tenplets (for example 10 ~ (a,b,Q), where 
a > b > Q is the charge of the first generation, b of 
the second, zero of the third) while no charge differ- 
ences appear in the 5 (e. g. 5 ~ (0, 0, 0)). The ob- 
served fact that the up-quark mass hierarchies are more 
pronounced than for down-quark and charged leptons is 
in agreement with this assignment. The S U{5) gener- 
ators act vertically inside one generation, whereas the 
U{l)fN charges are different horizontally from one gen- 
eration to the other. If, for a given interaction vertex, 
the U(1)fn charges do not add to zero, the vertex is for- 
bidden in the symmetric limit. However, the U{l)fN 
symmetry (that one can assume to be a gauge sym- 
metry) is spontaneously broken by the VEVs v/ of a 
number of flavon fields with non-vanishing charge and 
GUT-scale masses. Then a forbidden coupling is res- 
cued but is suppressed by powers of the small parame- 
ters A = Vf/M, with M a large mass, with the exponents 
larger for larger charge mismatch. Thus the charges fix 
the powers of A, hence the degree of suppression of all 
elements of mass matrices, while arbitrary coefficients 
kjj of order 1 in each entry of mass matrices are left un- 
specified (so that the number of parameters exceeds the 
number of observable quantities). A random selection 
of these kjj parameters leads to distributions of result- 
ing values for the measurable quantities. For anarchy 
(A) the mass matrices in the leptonic sector are totally 
random, while in the presence of unequal charges dif- 
ferent entries carry different powers of the order param- 
eter and thus some hierarchies are enforced. The em- 
bedding of Anarchy in the S U{5) ® U{1)fn context al- 
lows to implement a parallel treatment of quarks and 
leptons. Within this framework there are many variants 
of these models: fermion charges can all be nonnega- 
tive with only negatively charged flavons, or there can 
be fermion charges of different signs with either flavons 
of both charges or only flavons of one charge. In mod- 
els with no see-saw, the 5 charges completely fix the 
hierarchies (or anarchy, if the case) in the neutrino mass 
matrix. If Right-Handed (RH) neutrinos are added, they 
transform as S U(5) singlets and can in principle carry 
U{l)fN charges, which also must be all equal in the an- 
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archy case. With RH neutrinos the see-saw mechanism 
can take place and the resuhing phenomenology is mod- 
ified. In Ref. li28i] . given the new experimental results, 
we have made a reappraisal of Anarchy and its variants 
within the (SUSY) SU(5) x U(l)m GUT framework. 
Based on the most recent data we argue that the Anar- 
chy ansatz is probably oversimplified and, in any case, 
not compelling. In fact, suitable diff'erences of U(l)fff 
charges, if also introduced within pentaplets and sin- 
glets, lead to distributions that are in better agreement 
with the data with the same number of random param- 
eters as for Anarchy. The hierarchy of quark masses 
and mixing and of charged lepton masses in all cases 
impose a hierarchy defining parameter of the order of 
Ac = sin 6c, with 6c being the Cabibbo angle. The 
weak points of Anarchy (A) are that with this ansatz all 
mixing angles should be of the same order, so that the 
relative smallness of So ~ o{Ac) is not automatic. Sim- 
ilarly the smallness of r = Am^^^^^/Am^,,,, is not easily 
reproduced: with no See-Saw r is of o(l), while in the 
See-Saw version of Anarchy the problem is only par- 
tially alleviated by the spreading of the neutrino mass 
distributions that follows from the product of three ma- 
trix factors in the See-Saw formula. An advantage is 
already obtained if Anarchy is only restricted to the 23 
sector of leptons as in the A^^ model (in the notation of 
Ref. |i,28,l ). In this case, with or without See-Saw, So is 
naturally suppressed and, with a single fine tuning one 
gets both 012 large and r small (this model was also re- 
cently rediscussed in Ref. 1,291 ). Actually in Ref . 128 1 
we have shown that, in the no See-Saw case, a very good 
performance is observed in a new model, the H model, 
where Anarchy is also relaxed in the 23 sector In the 
H model, by taking a relatively large order parameter, 
one can reproduce the correct size for all mixing angles 
and mass ratios. Alternatively, in the See-Saw case, we 
have shown that the freedom of adopting RH neutrino 
charges of both signs, as in the f A^^ model, can be used 
to obtain a completely natural model where all small 
quantities are suppressed by the appropriate power of 
A. In this model a lopsided Dirac mass matrix is com- 
bined with a generic Majorana matrix to produce a neu- 
trino mass matrix where the 23 subdeterminant is sup- 
pressed and thus r is naturally small with unsuppressed 
623. In addition 612 is large while 613 is suppressed. We 
stress again that the number of random parameters is 
the same in all these models: one coefficient of o((l) for 
every matrix element. Moreover, with an appropriate 
choice of charges, it is not only possible to reproduce 
the charged fermion hierarchies and the quark mixing, 
but also the order of magnitude of all small observed 
parameters can be naturally guaranteed. In conclusion. 



we agree that models based on chance are still perfectly 
viable, but we consider Anarchy a particularly simple 
choice perhaps oversimplified and certainly not com- 
pelling and we have argued in favour of less chaotic so- 
lutions. 

Anarchy and its variants, all sharing the dominance 
of randomness in the lepton sector, are to be confronted 
with models with a richer dynamical structure, in partic- 
ular those based on discrete flavour groups (for a review, 
see, for example, Ref. f,30il ). After the measurement of 
a relatively large value for ^13 there has been an intense 
work to interpret these new results along different ap- 
proaches and ideas. An updated list of references can be 
found, for example, in the recent papers Refs. \ 51]. 
Among the models with a non trivial dynamical struc- 
ture those based on discrete flavour groups were mo- 
tivated by the fact that the data suggest some special 
mixing patterns as good first approximations like Tri- 
Bimaximal (TB) or Golden Ratio (GR) or Bi-Maximal 
(BM) mixing, for example. The corresponding mix- 
ing matrices all have sin^6'23 = 1/2, sin^^o = 0, 
values that are good approximations to the data (al- 
though less so since the most recent data), and differ 
by the value of the solar angle sin^6'i2. The observed 
sin^ 612, the best measured mixing angle, is very close, 
from below, to the so called Tri-Bimaximal (TB) value 
of sin^ 6*12 = 1 /3. Alternatively, it is 
also very close, from above, to the Golden Ratio (GR) 
value H Ell IMS sm^6n = ^ = 577! ~ 
where (p - (I + V5)/2 is the GR (for a different con- 
nection to the GR, see Refs. SEll]). On a diff'erent 
perspective, one has also considered models with Bi- 
Maximal (BM) mixing, where sin^fl^ - 1/2, i.e. also 
maximal, as the neutrino mixing matrix before diago- 
nalization of charged leptons. One can think of models 
where a suitable symmetry enforces BM mixing in the 
neutrino sector at leading order (LO) and the necessary, 
rather large, corrective terms to 612 arise from the diag- 
onaUzation of the charged lepton mass matrices (a list 
of references can be found in Ref. fsO*]). Thus, if one or 
the other of these coincidences is taken seriously, mod- 
els where TB or GR or BM mixing is naturally predicted 
provide a good first approximation (but these hints can- 
not all be relevant and it is well possible that none is). 
The corresponding mixing matrices have the form of ro- 
tations with fixed special angles. Thus one is naturally 
led to discrete flavour groups. Models based on dis- 
crete flavour symmetries, like A4 or 5 4, have been pro- 
posed in this context and widely studied. In these mod- 
els the starting Leading Order (LO) approximation is 
completely fixed (no chance), but the Next to LO (NLO) 
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corrections still introduce a number of undetermined pa- 
rameters, although in general much less numerous than 
for U{1)fn models. These models are therefore more 
predictive and typically, in each model, one obtains re- 
lations among the departures of the three mixing angles 
from the LO patterns, restrictions on the CP violation 
phase 6cp, mass sum rules among the neutrino mass 
eigenvalues, definite ranges for the neutrinoless beta de- 
cay effective Majorana mass and so on. 

In the following we will mainly refer to TB or BM 
mixing which are the most studied first approximations 
to the data. A simplest symmetry that, at LO, leads to 
TB is A4 while BM can be obtained from 5 4. Starting 
with the ground breaking paper in Ref. (a^, A4 mod- 
els have been widely studied (for a recent review and a 
list of references, see Ref. 14911 ). At LO the typical A4 
model (like, for example, the one discussed in in Ref. 
14311) leads to exact TB mixing. The LO approxima- 
tion is then corrected by non-leading effects. Given the 
set of flavour symmetries and having specified the field 
content, the non-leading corrections to TB mixing, aris- 
ing from higher dimensional effective operators, can be 
evaluated in a well-defined expansion. In the absence 
of specific dynamical tricks, in a generic model all three 
mixing angles receive corrections of the same order of 
magnitude. Since the experimentally allowed depar- 
tures of O12 from the TB value, sin^ O12 = 1/3, are small, 
numerically not larger than 0{A^) where Ac = sin 6c, it 
follows that both 613 and the deviation of 623 from the 
maximal value are also expected to be typically of the 
same general size. The same qualitative conclusion also 
applies to A5 models with GR mixing. This generic pre- 
diction of a small ^13, numerically of 0{A^), can now be 
confronted with the present precise experimental value. 
The central value sin 613 ~ 0.15, from Table 2, is be- 
tween 0(Al) ~ (9(0.05) and 0(Ac) ~ 0(0.23). Since Ac 
is not that small, this gap is not too large and one can 
argue that models based on TB (or GR) mixing are still 
viable with preference for the lower side of the experi- 
mental range. 

Of course, one can introduce some additional theoret- 
ical input to improve the value of 613. In the case of A4, 
one particularly interesting example is provided by the 
Lin model (s^] (see also Ref. [EllsTl), formulated be- 
fore the T2K, MINOS, DOUBLE CHOOZ, Daya Bay 
and RENO results. In the Lin model the A4 symme- 
try breaking is arranged, by suitable additional Z„ par- 
ities, in a way that the corrections to the charged lep- 
ton and the neutrino sectors are kept separated not only 
at LO but also at next-to-leading order (NLO). As a 
consequence, in a natural way the contribution to neu- 
trino mixing from the diagonalization of the charged 



leptons can be of O(A^), while those in the neutrino 
sector of 0(Ac). In addition, in the Lin model these 
large corrections do not affect 612 and satisfy the relation 
sin^ 023 -1/2+ 1/ VZcos^cpI sin0i3|, with 6cp being 
the CKM-like CP violating phase of the lepton sector. 
Thus, in the Lin model the NLO corrections to the solar 
angle 612 and to the reactor angle 813 are not necessarily 
related. Note that, for 623 in the first octant, the sign of 
cos 6cp must be negative. 

Alternatively, one can think of models where, be- 
cause of a suitable symmetry, BM mixing holds in the 
neutrino sector at LO and the corrective terms for 6*12, 
which in this case are required to be large, arise from 
the diagonalization of charged lepton masses. These 
terms from the charged lepton sector, numerically of or- 
der 0{Ac), would then generically also affect 613 and the 
resulting angle could well be compatible with the mea- 
sured value. An explicit model of this t ype based on 
the group 5 4 has been developed in Ref. ll45ll (see also 
Refs. 11461 147I ksll ). An important feature of this model 
is that only 612 and 613 are corrected by terms of 0(Ac) 
while 023 is unchanged at this order. This model is com- 
patible with present data and clearly prefers the upper 
range of the present experimental result for 613. 

In Ref. 15511 we discuss three possible classes of mod- 
els: 1) typical A4 models where 813 is generically ex- 
pected to be small, of the order of the observed depar- 
tures of 012 from the TB value, and thus with preference 
for the lower end of the allowed experimental range. 
2) special A4 models, like the Lin model, where 0^3 is 
made independent of the deviation of sin^ 0\2 from the 
TB value 1/3 and can be as large as the upper end of 
the allowed experimental range. In the same paper we 
discuss a general characterization of these special A4 
models where the dominant corrections to TB mixing 
do not arise from the charged lepton sector but from the 
neutrino sector. 3) Models where BM mixing holds in 
the neutrino sector and large corrections to 6*12 and ^13 
arise from the diagonalization of charged leptons. The 
value of 013 is naturally close to the present experimen- 
tal range. In each of these possible models the dominant 
corrections to the LO mixing pattern involve a number 
of parameters of the same order of magnitude, ^. We 
discuss the success rate corresponding to the optimal 
value of ^ for each model, obtained by scanning the pa- 
rameter space according to a similar procedure for all 
three cases. We argue that, while the absolute values of 
the success rates depend on the scanning assumptions, 
their relative values in the three classes of models, pro- 
vide a reliable criterium for comparison. We find that, 
for reproducing the mixing angles, the Lin type models 
have the best performance, as expected, followed by the 
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typical A4 models while the BM mixing models lead to 
an inferior score, as they can well reproduce the size of 
6*13 but most often fail to reproduce the correct value of 
9i2- We also discuss the conditions for cos^cp peaking 
around - 1 . 



In Ref. issi we have also discussed the implica- 
tions for lepton flavour-violating (LFV) processes of 
the above three classes of possibilities, assuming a su- 
persymmetric context, with or without See-Saw. The 
present bounds on LFV reactions pose severe con- 
straints on the parameter space of the models (for a 
recent general analysis on model-independent flavour 
violating eff'ects in the context of flavour models, see 
Ref. 15611 ). In particular, we refer to the recent im- 
proved MEG result 115 711 on the ^ ^ ey branching ra- 
tio, Brin ey) < 2.4 x 10"'^ 95% C.L. and to 
other similar processes like t — > (e or ^)y. One expects 
that lepton flavour-violating processes may also have a 
large discriminating power in assessing the relative mer- 
its of the diff'erent models. In Ref. ll55ll we have studied 
this issue by adopting the simple CMSSM framework. 
While this very constrained version of supersymmetry 
is rather marginal after the results of the LHC searches, 
more so given that the Higgs mass is around mn = 125 
GeV, we still believe it can be used for indicative pur- 
poses as in this case. We find that the most constrained 
versions are the models with BM mixing at LO where 
relatively large corrections directly appear in the ofl"- 
diagonal terms of the charged lepton mass matrix. The 
A4 models turn out to be the best suited to satisfy the 
experimental bounds, as the non-diagonal charged lep- 
ton matrix elements needed to reproduce the mixing an- 
gles are quite smaller. An intermediate score is achieved 
by the models of the Lin type, where the main correc- 
tions to the mixing angles arise from the neutrino sec- 
tor and the non-diagonal charged lepton matrix elements 
are smaller. Overall, among the discrete flavour group 
models that we have studied, the A4 models emerge well 
from the analysis of Ref. fss'land in particular those of 
the Lin type lead to a natural description of the data. 
As for the regions of the CMSSM parameter space that 
are indicated by our analysis, the preference is for small 
tan p and large SUSY masses (at least one out of otq and 
mi/2 must be above 1 TeV). Given that SUSY has not 
been found at the LHC the preference for heavy SUSY 
masses does not pose a problem to these flavour models, 
except that, as a consequence, it appears impossible, at 
least within the CMSSM rigid framework, to satisfy the 
MEG bound and simultaneously to reproduce the muon 
g - 2 isi discrepancy iHliolleil]. 
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